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The connection between the asymptotic states of the S'-matrix and the energy spectrum in a finite volume 
is investigated in the context of a ANn interaction model that can be solved exactly. In particular, Luscher's 
method, which proposes a definite relationship between the elastic scattering phase shift and the finite- volume 
energy levels, is shown to exhibit a deviation from the correct resonance position that behaves approximately 
like 1/L in the range of volumes of interest. A new method is proposed, whereby the energy levels from a 
lattice calculation are used to determine the parameters of the model. The phase shift may then be evaluated 
using an infinite- volume t-matrix. The model-dependence of the method is investigated. At box sizes that are 
currently considered quite large, model-dependent effects are estimated to be an order of magnitude smaller 
than the deviations associated with Luscher's method. The extraction of the resonance position is also shown to 
be relatively independent of L. 
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Recent advances in lattice QCD have demonstrated the abil- 
ity to extract the low-lying tower of excitation energies re- 
liably in baryon correlators [1-3]. Nevertheless, the identi- 
fication of resonance parameters (e.g. masses, widths, etc.) 
from multi-hadron interactions remains an ongoing challenge 
of current and active research [4—14]. In particular, the cal- 
culation of the discrete energy levels of a hadronic system on 
a finite-volume lattice needs to be related to the continuous 
spectra measured by scattering experiments. One approach, 
Luscher's method [15], allows such a relation to be made by 
geometrically matching the asymptotic behaviour of the S'- 
matrix onto the energy spectrum of the toroidal topology of 
the lattice. The standard implementation of Luscher's method 
neglects corrections to energy levels that are exponentially 
suppressed in the lattice volume. In this Letter, an exactly 
solvable model is used to investigate the role of these terms in 
the extraction of resonance parameters. 

Luscher's method entails a fixed relationship between the 
scattering phase shift S, and the rth energy level in a finite 
volume (in keeping with Luscher's notation, r is an integer): 



r 7r — S(k r ) = 4>(q r ). 
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q r = k r L/ (2ir) are the (dimensionless) lattice momenta, and 
4>(q) is a geometric function related to the three-dimensional 
Zeta function, defined by: 
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A simple minimal subtraction scheme has been chosen to reg- 
ulate the divergence in the Zeta function. 

Luscher's method will be tested in the context of a well- 
studied model describing the decay of the A-baryon into A^7r; 
an archetypal example of a baryon resonance. It will be shown 
that Luscher's method entails a deviation from the infinite- 
volume resonance position of order ~ 1/L, for moderate and 



even large volumes. A new method is proposed, in which 
energy levels from the lattice are fit to those of the model, 
and the phase shift then calculated using the infinite-volume 
t-matrix. 

In the following finite-volume heavy-baryon model, a ma- 
trix Hamiltonian, H = Hq + Hj, is defined to include the 
interaction between the A-baryon and the pion-nucleon sys- 
tem. The rows and columns of H represent the momentum 
values available to the pion (fci = 2ir/L, . . .). In this model, 
the Nir states are only coupled to the A-baryon. The free 
Hamiltonian, Hq, takes the form: 
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where lo v (k n ) = -\/fc 2 + m 2 , and the momenta at finite vol- 
ume satisfies the condition: fc 2 = (^7) 2 n, for a squared in- 
teger n = n 2 + n 2 + n 2 . The matrix H has non-zero diag- 
onal entries representing the non-interacting energy of a pion 
at certain values of momentum, and the nucleon recoil van- 
ishes in the heavy-baryon limit. Note that the cubic symmetry 
of the system means that the dimension of H is significantly 
smaller; each non-zero n occurs only once. The first element 
of the Ho includes a bare resonance mass relative to the nu- 
cleon mass, An, which is chosen so that the phase shift passes 
through 90° at 292 MeV in the infinite-volume limit. 

The interacting Hamiltonian, Hj, contains couplings, gAN, 
in the top row and left-most column: 
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The A -H> Ntt coupling, g% n N {k n ), is defined in terms of 
known phenomenological parameters f v = 92.4 MeV, and 
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C = —1.52, as derived from the SU(6) flavour-symmetry re- 
lation, and must also contain a normalization factor yj 63(71), 
representing the number of ways of summing the squares of 
three integers to equal n: 
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and u(k n ) = 1 
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Note that a regulator, u(k), has been introduced into the 
coupling term. This will allow the form of the energy levels 
to be related to the regularized one-pion loop integral occur- 
ring in the t-matrix, as shown below. (For detailed discussions 
of finite-range regularization, see Refs. [16-20]). For the pur- 
poses of this model, a dipole regulator with A = 0.8 GeV is 
chosen. It will be shown, however, that the form of the reg- 
ulator has little impact on the consistency of the final result: 
extracting the correct resonance position. 

One extension of this model would be to include a direct 
Ntt «-» Ntt interaction, which could be viewed as an approx- 
imation to the famous Chew-Low interaction. For example, 
one could add additional terms in the interaction matrix Hj : 



9CL(ki,kj) = -gcL(ki)g C L(kj) 
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In the Cloudy Bag Model [21, 22], The couplings gcL and 
are related by the equation: gcL(ki) = ng A n N (ki), with 
a ratio, n = -\/25/18, estimated from SU(6) symmetry. In the 
case of the t-matrix, the Lippmann-Sch winger equation with a 
rank-2 separable potential can be solved as per Mongan [23]. 
Suffice to say, this additional interaction will be considered 
only in the final analysis, as a test of model-independence. 

Firstly, consider only the A f-> Ntt interaction. The eigen- 
values of the matrix H correspond to the different energy lev- 
els in a finite volume, and may be calculated for a range of 
L values. The ten lowest energy levels as a function of L are 
shown in Fig. 1, with the non-interacting energies shown as 
dotted lines. It will be useful to have a reliable method for 
extracting phase shifts in this region, even if larger box sizes 
were feasible in lattice QCD. 

In the case of this model, the phase shift can, in fact, be 
calculated exactly, thus providing a suitable test for the ability 
Liischer's method to reproduce the correct phase shift. The 
solution to the eigenvalue problem for H takes the following 
form (provided uo n {k n ) ^ A): 
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which mirrors the i-matrix calculation for elastic AT7r- 
scattering with a A-baryon intermediate. The relationship 
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FIG. 1. (color online). The ten lowest energy levels from the A7V7T model 
described in Eqs. (4)— (9) (solid lines), and the corresponding non-interacting 
energies (dotted lines). 



between the phase shift and the on-shell quantity, T 

t(k,k;E+), is: 
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The bare resonance mass, A , may be tuned so that the final 
resonance energy matches the physical AN mass-splitting: 
292 MeV, and consequently fed into Eq. (4). 

The A -H> Ntt one-pion loop integral in a finite volume 
is equal to the sum in Eq. (11), matching the effective field 
theory calculation in the weak coupling limit, Ea —} Aq: 
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At infinite volume, the real part of the loop integral takes the 
form: 
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where V indicates that a principal value integral must be per- 
formed. By solving Eq. (12), the phase shift, S, may be plot- 
ted as a function of energy, E — ui n (k), to obtain the Breit- 
Wigner-like curve shown in Fig. 2. 

In order to test Liischer's formula for the AA^7r model, the 
phase shift is calculated from the eigenvalues obtained from 
H using Eq. (1). This phase shift can be compared to the 
infinite-volume phase shift calculated from Eq. (12). A plot 
of these estimates, 5(k r ), for a range of volumes is shown 
alongside the infinite-volume phase shift as a function of E in 
Fig. 3. It is clear that the convergence to the infinite-volume 
phase shift as a function of L is quite slow, with an error of 
approximately 20 MeV in the resonance position, £V 0S , on a 5 
fm box (using a linear interpolation between the two energies 
closest to E les ). 
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FIG. 2. (color online). The Breit-Wigner-like phase shift S associated with 
elastic 7V7r-scattering via a A-baryon intermediate state, plotted against ex- 
ternal energy E, as calculated from the on-shell i-matrix of Eq. (12). 
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FIG. 3. (color online). Estimates of the phase shift obtained from Liischer's 
formula applied to the ANtt model for a range of L values. The dot-dashed 
lines that connect the discrete points are merely to guide the eye. The infinite- 
volume phase shift is also shown as a solid line. 



By applying Liischer's formula from Eq. (1), the corre- 
sponding estimate of the energy levels may be calculated. 
The energy difference between the exact energy levels from 
the model and the levels obtained from Liischer's formula is 
shown for the five lowest-lying states in Fig. 4. The energy 
difference is exponentially suppressed for large values of L. 

The behaviour of E ICS as a function of L is of great interest. 
This provides an indication of the robustness of the method 
for extracting an estimate of the true phase shift. The result of 
obtaining E Tea , using a linear and also a Breit-Wigner inter- 
polation, is shown as a function of 1/L in Fig. 5. Firstly, there 
are significant oscillations in the behaviour of E res even up to 
box sizes of 7 fm. This is a consequence of the behaviour of 
the geometric <f> function. Furthermore, the large L behaviour 
is linear in 1/L; though it is reasonable to expect that expo- 
nential suppression will dominate for larger box sizes than 40 
fm. A linear extrapolation of the resonance position in 1/L 
is possible, but Fig. 5 indicates that a reliable result requires 
a box size of L > 10 fm. Clearly, these box sizes far exceed 
feasibility for a current lattice QCD calculation. For presently 
conceivable lattice sizes (L < 5 fm), the error in the resonance 
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FIG. 4. (color online). The energy differences 8E between the exact levels 
from the ANtt model and the levels obtained from Liischer's formula using 
the exact phase shifts from the t-matrix. Energy level 1 corresponds to the 
lowest lying state. 
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FIG. 5. (color online). The resonance energy E TCB plotted against 1/L, 
using a linear and a Breit-Wigner interpolation of the discrete phase shift 
obtained using Eq. (1). The experimental value is marked with a square. The 
results from fitting pseudodata with a different regulator, a Gaussian with A = 
0.5, 0.6 & 0.8 GeV, are plotted. The results of fitting pseudodata containing 
anon-resonant background interaction (see Eq. (10)) are also shown, treating 
the Gaussian scale A as a fit parameter. 



position obtained through either a linear or a Breit-Wigner in- 
terpolation to the point where 5 = 90°, yields an error of order 
10%. For this reason, we have been led to consider an alter- 
native method for extracting an estimate of the phase shift. 

Motivated by the general result that the i-matrix near any 
resonance is separable [24], the idea of fitting the energy lev- 
els from a lattice calculation to a separable model is explored. 
By treating the coupling constant, \A, and the bare resonance 
mass, Ao, as fit parameters, the energy levels between the lat- 
tice calculation and the previously established ANtt model 
can be matched. A new estimate of the exact phase shift can 
be then calculated from the t-matrix calculation of Eq. (12) 
using the parameters of the fit. 

In order to test this idea, the energy levels for the ANtt 
model are generated using a dipole regulator with A = 0.8 
GeV, and a value of Ao tuned to reproduce the experimental 
resonance position. These energy levels are treated as pseu- 
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dodata. A model is then chosen using a different regulator 
function, such as a Gaussian regulator. By minimizing the x 2 
function between the two sets of energy levels and obtaining 
fit values for xa and A , exact infinite-volume phase shifts 
may be calculated for a variety of values of L. Only energy 
levels in the region where the model is valid, near the reso- 
nance (E r < 500 MeV), are used when matching two sets of 
energy levels. The resonance position obtained using this pro- 
cedure is shown in Fig. 5 as a function of 1 /L, for a variety 
of suitable values of the Gaussian regulator mass. Clearly, re- 
fitting the pseudodata with a different regulator function still 
leads to a result that is stable with respect to L. Furthermore, 
the extracted resonance position is within a few MeV of the 
actual resonance position. This is encouraging, because it 
suggests that approximating the underlying physics of a lat- 
tice calculation with a regulator, sensibly chosen, will lead to 
a result that is not highly dependent on the features of the par- 
ticular form of the regulator function. 

One might argue that the simple model including only 
A o Nn couplings is ideally suited to the reconstruction 
method proposed here. For that reason, the following case 
was considered, in which pseudodata were generated from a 
Hamiltonian model that now includes the effective contact in- 
teraction from Eq. (10), again using a dipole regulator with 
A = 0.8 GeV. The energy levels were fit to the more basic 
model, with only the A <H> Nit coupling, using a Gaussian 
regulator, and treating A as an additional fit parameter. In this 



case, the fitted value of the Gaussian regulator is Afi t = 0.42 
GeV, leading to a stable extraction of the resonance energy, 
as shown in Fig. 5. This is remarkable, as it shows that the 
bare resonance mass, the coefficient of the interaction and the 
finite-range scale A are collectively able to mock-up the miss- 
ing Chew-Low interaction. That is, more complicated inter- 
actions are incorporated into the fit degrees of freedom. 

In summary, it has been demonstrated that care must be 
taken in matching the asymptotic states of elastic scattering 
to the energy spectrum in a finite volume. The finite range 
of the hadronic interactions entails a ~ 1/L deviation in the 
extracted resonance position using Liischer's method. A new 
method is proposed in which a A Ntt interaction model is con- 
structed, allowing the bare resonance position and the cou- 
pling strength to be free parameters. By matching the energy 
levels of an ab initio finite-volume calculation, such as lat- 
tice QCD, with those of the model, the fitted parameters can 
be used to reconstruct the exact phase shift using the infinite- 
volume t-matrix calculation. The results show that the reso- 
nance position of the new phase shift is both accurate, and rel- 
atively independent of L across a wide range of box sizes, in- 
cluding those used in contemporary lattice QCD simulations. 
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